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Abstract 

The  side-by-side  computational  and  experimental  combustion  research  has  made  remarkable  progress  in 
unifying  concepts  of  theoretical  foundation  for  canonical  laminar  flame  phenomena.  It  is  logical  to  adopt  the 
same  approach  to  turbulent  combustion  by  explore  this  cutting  edge  technology  [1-3].  However  the  turbulent 
structure  in  combustion  presents  a  formidable  challenge.  Turbulence  is  the  most  complex  fluid  physics 
phenomenon  of  nonlinear  dynamics  associated  with  multi-scales  in  time  and  space.  The  drastic  difference  in 
flow-field  structures  between  the  laminar  and  turbulent  counter- flowing  fuel  injection  is  clearly  illustrated  in 
figure  1.  As  a  consequence,  it  has  long  been  recognized  as  the  pacing  item  in  computational  fluid  dynamics 
(CFD)  [1],  This  urgent  scientific  need  for  technologic  advancement  is  far  from  being  met.  Two  fundamental 
approaches  have  been  developed:  the  Direct  Numerical  Simulation  (DNS)  [4,5]  and  the  Large  Eddy  Simulation 
(LES)  techniques  [6,7],  According  to  the  Kolmogorov  micro-scales 
formulation,  the  ratio  of  the  extreme  length  scales  is  proportional  to 
three-quarter  power  of  the  characteristic  Reynolds  number,  Re3  4  of 
a  turbulent  flowfield.  To  resolve  all  scales  for  three-dimensional 
motion,  the  number  of  grid  points  required  becomes  astronomical  to 
be  a  practical  modeling  and  simulation  tool  [1,4].  In  principle,  LES 
can  accurately  compute  the  contribution  of  the  large  energy¬ 
carrying  structures  for  momentum  and  energy  transfer  if  numerical 
error  is  minimized.  Meanwhile,  the  fine-scale  structures  of 
combustion  depend  upon  intermolecular  interactions  such  as 
transport  properties  and  chemical  reaction  beneath  the  molecular 
scales,  which  must  be  modeled  [2-11].  Thus,  the  required 
computational  resource  is  also  prohibitive  for  practical  application 
at  the  present  time. 

An  alternative  description  for  turbulent  flames  is  the  concept  of  laminar  flamelet,  which  has  been  introduced  by 
Peters  [2]  as  an  asymptotic  concept.  This  idea  is  built  on  the  premises  that  combustion  takes  place  within 
asymptotically  thin  layers  embedded  in  the  turbulent  flow.  Therefore  the  local  instantaneous  composition  of 
reactants  and  temperatures  of  the  non-premixed  process  could  be  identically  modeled  as  that  in  a  stretched 
laminar  flame.  A  necessary  condition  for  the  flamelet  concept  to  be  validated  is  that  the  characteristic 
dimension  of  the  reaction  zone  must  be  thinner  than  the  Kolmogorov  length  scale. 

However,  Bilger  has  shown  the  reaction  zone  of  the  flamelet  is  not  necessary  asymptotically  thin  due  the 
reverse  reactions  [3],  It  is  known  from  numerous  calculations;  the  minimum  reaction  zone  for  hydrogen  flames 
has  a  thickness  from  0.5  to  1.0  mm  that  is  aerodynamically  stretched  [8-10].  Recent  numerical  simulations  of 


vortex-flame  interaction  by  Katta  el  al  [8]  indicate  that  when  the  vortices  are  smaller  than  millimeter  in 
dimension  will  not  stretch  a  laminar  non-premixed  flame  but  destroy  the  flame  structure  in  the  process.  This 
finding  is  very  interesting  it  also  raises  the  question  whether  the  laminar-flamelet  theory  is  applicable  to 
turbulent  combustion  [9],  In  order  to  further  substantiate  the  discovery,  a  local  highly  refined  grid  spacing 
algorithm  of  a  spectral-like  numerical  resolution  together  with  a  refined  diffusion  model  becomes  necessary  [8- 
10].  Since  the  chemical  reaction  occurs  beneath  the  molecular  scales  which  is  irresolvable  by  the  classic  theory 
of  gas  kinetics  and  must  be  augmented  by  models  of  chemical  kinetics.  Fortunately,  this  scientific  disciplinary 
has  been  developed  by  combustion  communities,  the  reaction  rate  of  the  chemical  kinetics  and  the  molecular 
and  atomic  diffusion  model  can  be  accessed  through  the  Cyber  space  [1 1,12].  The  present  research  is  therefore 
required  only  to  focus  on  the  development  of  a  high  resolution  algorithm  that  is  supported  by  extremely  fine 
local  mesh  spacing  for  combustion  simulation. 

The  formulation  of  Gauss  quadrature  is  known  for  its  high  numerical  resolution  in  performing  integration 
between  intervals  of  a  real  function.  The  analytic  formation  is  expressed  by  the  Newton’s  division-difference 
form  on  unequal  roots  of  the  adopted  orthogonal  polynomial  [13,14],  The  basic  formulation  however  is 
applicable  to  derivative  calculations  for  any  function  within  the  locally  refined  domain  [15-19].  The  calculated 
derivative  by  this  approach  can  yield  a  spectral-like  resolution  more  accurate  than  most  high-order  numerical 
schemes  and  the  accuracy  is  controlled  solely  by  the  degrees  of  the  adopted  polynomial.  For  this  reason,  the 
ability  to  refine  local  grid  spacing  within  an  investigated  domain  using  an  orthogonal  polynomial  is  also 
unsuppressed  [18,19].  In  general,  the  solutions  between  immediately  adjacent  domains  are  discontinuous.  In 
applying  this  approach  to  solve  hyperbolic  partial  differential  equations,  the  conservation  law  at  the  interface 
boundary  between  refined  domains  requires  only  to  balance  the  outward  normal  component  of  the  flux  vector. 
The  solution  is  permitted  to  be  discontinuous  at  the  domain  interface  to  satisfy  the  conservative  property,  and 
the  approximate  Riemann  approach  is  well-posed  [20,21],  However  in  solving  the  diffusion  dominant  partial 
differential  equations  for  combustion,  the  function  continuity  and  adjoint  consistency  of  solutions  between 
refined  domains  must  be  maintained  across  the  interface  and  this  issue  has  been  known  at  the  very  beginning  of 
the  development  of  the  polynomial  refinement  algorithms  [15].  A  wide  range  of  procedures  have  been 
investigated  from  the  penalty  method  [22]  to  the  more  recent  studies  on  the  recovery  polynomial  basis 
approaches  [18,19,23,24],  The  essential  idea  is  rest  on  an  effective  coupling  of  the  discontinuous  solutions 
across  the  interface  between  locally  refined  domains  which  permits  a  very  large  class  of  solving  schemes  to  be 
utilized  [18,19].  This  issue  has  been  one  of  the  focused  areas  of  the  present  intensive  basic  research. 

The  inherent  deficiency  of  the  local  polynomial  refinement  using  Gauss  Quadrature  formulation  at  the  interface 
of  locally  refined  domains  is  remedied  by  a  differential  reconstruction  scheme  [18,19],  The  continuous  solution 
at  the  interface  for  combustion  simulation  is  enforced  by  this  approach  for  the  diffusion  dominant  equation.  A 
multiple-point,  unequal  spaced,  differencing  reconstruction  approximation  in  preserving  the  spectral-like 
accurate  and  to  increase  computationally  efficient  has  been  accomplished.  It  is  found  that  the  shared  first  and 
second  derivatives  at  the  interface  can  be  computed  from  the  roots  of  polynomials  of  adjacent  domains  and 
effectively  retains  the  numerical  resolution  in  comparison  with  the  recovery-based,  L2  norm  projection,  penalty, 
and  reconstruction  methods  in  applications  [18,19,24],  The  computed  flame  thickness  agrees  well  with  results 
using  a  systematic  mesh  refinement  approach  and  experimental  observation  [8-10,18,19], 


The  intrinsic  behavior  of  the  present  algorithm  in  practical 
applications  has  been  successfully  demonstrated  for  the 
counter-flowing,  air-hydrogen  combustion  as  shown  in  Figure 
2.  The  highly  stretched  flame  that  approaching  the  extinct 
state  is  successfully  simulated;  the  thin  flame  structure  is 
accurately  captured  by  the  developed  spectral  method  through 
the  diminishing  residual  formulation  [1],  Meanwhile, 
improvement  has  also  been  made  to  increase  the  convergence 
rate  by  addressed  the  flux  vector  balancing  and  continuity 
preserving  condition  across  the  interfaces  of  the  immediately 
adjacent  computational  domains  [10,18,19], 
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Figure  2.  Detailed  structure  of  stretched  counter¬ 
flowing  air-hydrogen  flame 


In  order  to  cover  a  very  wide  range  of  computational  simulations  from  the  diffusive  counter-flowing  hydrogen 
flame  to  the  ethylene/JP7  scramjet  combustion,  a  system  of  robust  and  accurate  numerical  procedures  becomes 
a  critical  element  for  modeling  and  simulation.  The  overall  resolution  also  becomes  an  open  issue  when 
Riemann  problem  is  invoked  for  shock  capturing.  The  incurred  Gibbs  phenomenon  in  numerical  solutions  must 
be  controlled  by  either  adding  artificial  dissipation  or  using  a  post-process  with  low-pass  filters  [25],  The 
resultant  numerical  resolution  is  thereby  uncertain.  This  concerned  is  also  existed  when  applying  a  high  order 
WENO  scheme  to  solve  hyperbolic  conservation  laws  for  flowfield  with  shock  jumps,  the  improved  numerical 
behavior  is  not  derived  from  its  superior  accuracy  at  the  critical  points  of  the  solution,  but  it  rather  gains  through 
the  larger  weights  assigning  to  stencils  with  discontinuities  [26],  To  address  this  concern,  a  comparative  study 
of  the  adaptive  polynomial  local  refinement  and  the  compact-differencing  up  to  sixth-order  accuracy  has  been 
performed  to  assess  the  ability  in  defining  the  fine-structure  of  a  diffusive  flame  and  in  sustaining  supersonic 
combustion.  For  the  compact-differencing  scheme,  the  staged  one-order-lower  approximation  for  the  transition 
operator  between  the  interior  domain  and  the  far-field  boundaries  is  used  and  this  formulation  is  based  on  the 
summation-by-parts  energy  mode  analysis  [18,21,27]. 

The  comparative  investigation  of  the  local  grid  refinement  and  the  compact-differencing  schemes  are 
implemented  by  the  diminishing  residual  algorithm  [1,10,18,19].  In  this  procedure,  the  discretized  partial 
differential  system  is  first  cast  in  the  delta  form  in  time,  and  then  the  equations  are  further  split  into  different 
left-hand-side  (LHS)  and  right-hand-side  (RHS)  differencing  operator.  The  complete  governing  equation  at  any 
instance  is  retained  on  the  RHS  of  the  equation  and  is  calculated  by  the  high  resolution  numerical  schemes.  The 
LHS  is  required  only  to  maintain  computational  stability  to  drive  the  iterative  process  to  a  vanishing  temporal 
residual.  Whence  the  residual  has  reached  a  designated  criterion,  the  solution  will  yield  the  desired  spectral-like 
accuracy  correspond  to  the  implemented  local  polynomial  refinement  or  compact-difference  schemes.  A 
systematic  evaluation  of  the  high-resolution  schemes  in  the  high  local  gradient  domain  and  its  transient 
characteristic  to  the  far  field  has  been  successfully  accomplished  [18,19,28]. 

The  quantification  of  comparisons  is  first  given  by  the 
temperature  distribution  along  the  center  line  of  the  heat 
source  from  combustion  in  Figure  3.  The  levels  of  the 
resolved  temperature  are  affected  by  the  discontinuous  heat 
source  distribution;  the  temperature  by  the  windward  scheme 
show  numerical  oscillation  at  the  jumps  and  results  in  a 
much  lower  temperature,  980.73  K,  and  the  numerical 
oscillations  exhibit  a  magnitude  of  ±27.50  K  at  both  terminal 
locations.  There  also  is  a  difference  between  the  predictions 
by  the  polynomial  local  refinement  and  the  fourth-order 
compact-difference  schemes.  The  former  yields  a  value  of 
1120.00  versus  the  latter  of  1100.00  K.  At  the  edge  of 
discontinuous  heat  source  from  combustion,  the  jumps  are 
captured  by  the  polynomial  local  refinement  method  in  six 
unequal-spacing  root  points.  The  resolution  by  the  compact- 
difference  scheme  requires  two  more  equal-space  points. 

The  comparisons  of  the  first  and  the  second  spatial  derivatives  along  the  centerline  of  the  heat  source  from 
combustion  are  also  obtained  and  depicted  in  Figures  4  and  5.  In  spite  of  the  difference  in  the  supporting  mesh- 
point  number  density,  the  first  and  second  derivatives  by  the  polynomial  local  refinement  and  the  compact- 
difference  schemes  in  the  middle  of  the  global  grid  blocks  are  nearly  identical  to  each  other,  to  yield  the 
absolute  values  of  0.2780xl06  versus  0.2746xl06  and  0.3292xl09  versus  0.3135xl09  respectively.  The 
conventional  windward  approximation  underpredicts  the  derivation  by  23.6%  with  respect  to  the  high- 
resolution  procedures  also  shows  numerical  oscillations  approaching  and  leaving  the  abruptly  changing  heat 
source  locations.  In  short,  both  the  high-resolution  schemes  exhibit  a  desired  improvement  over  the  windward 
approximation  for  determining  the  maximum  attainable  temperature  along  the  heat  source  from  combustion  and 


K 

Figure  3.  Comparison  of  nuemrical  resolution  at  the 
heat  source  of  combustion 


suppressing  the  oscillatory  numerical  behavior.  These  improvements  in  numerical  resolution  are  essential  for 
investigating  the  detailed  flame  structure  with  fine-scale  vortical  interactions. 


X 

Fig.  4  Comparison  of  calculated  first  derivative  (ST/dx)  for 
the  two-dimensional  model  equation 


Figure  6  presents  the  final  assessment  for  the  developed  local  polynomial  refinement  procedure  based  on  the 
Gauss-Lobatto  formulation  by  comparing  the  hydroxide  concentration  in  molar  fraction  with  the  benchmark 
results  of  AFRL/RZ.  This  validating  base  is  generated  by  the  UNICORN  program  on  the  finest  mesh  available 
and  the  local  polynomials  refinement  is  achieved  by  16  global  grid  blocks  with  three  different  degrees  of 
polynomial  of  11,  25,  and  51.  The  highest  degree  polynomial  (51)  is  applied  only  in  the  8th  grid  block  (17.5  <  r 
<  20  mm).  The  immediately  adjacent  grid  blocks  to  the  highest  gradient  domain  are  assigned  the  25-degree 
polynomials  and  the  11 -degree  polynomials  farther  away  from  the  high  gradient  region.  The  flame  thickness 
determined  from  the  radical  component,  OH,  of  the  hydrogen-air  combustion  is  slightly  thinner  than  that  of  the 
molecular  component  of  hydrogen.  In  essence,  the  flame  thickness  varies  from  chemical  species  to  species 
according  to  the  associated  diffusion  velocity  and  reaction  rate.  In  essence,  the  flame  thickness  is  best 
characterized  by  the  thermal  thickness  of  the  flame.  From  the  computational  simulation,  the  counter- flow  jet 
velocity  exhibits  a  complex  pattern  which  decelerates  toward  the  flame  and  accelerates  at  the  outer  edge  of  the 
flame  then  finally  reaches  a  free-stream  stagnation  point  within  the  fame  of  the  counter-flow  combustion 
process.  The  complex  flame  structure  is  clearly  governed  by  the  molecular  diffusion  and  the  rate  of  species 
generation  and  recombination.  An  in-depth  analysis  can  be  achieved  only  by  examining  all  chemical  species 
within  the  flame. 


The  agreement  between  solutions  of  the  present 
computational  simulations  with  AFRL  benchmark  is 
excellent.  The  maximum  discrepancy  between  the  two 
results  is  4. 172x1  O'3  in  molar  fraction  of  hydroxide,  and 
is  less  than  0.23xl0'2  mm  in  predicted  flame  thickness  of 
2.07  mm.  However,  the  finest  grid  spacing  for  the 
locally  refined  computational  domain  is  7. 194x1  O'5  mm 
and  the  detailed  flame  structure  is  described  by  29 
polynomial  roots.  It  is  important  to  note  that  the  solution 
using  Gauss-Lobatto  formula  can  still  be  refined  through 
the  adaptation  of  an  even  higher  and  unlimited  degree  of 
polynomial.  For  the  present  simulation  condition,  a  total 
of  243  root  points  of  the  Legendre  polynomial  are  used. 
The  validating  agreement  between  results  is  anticipated 
because  of  the  shared  kinematic  and  thermodynamic 
data  from  the  UNICORN  program.  The  most  important 
and  meaningful  conclusion  can  be  drawn  from  this  comparison  is  that  the  Gauss-Lobatto  formulation  has  the 
desired  resolution  characteristic  for  the  intermediate  wave  number  within  a  flame.  It  offers  a  potentially  great 
advantage  in  practical  application  for  combustion  simulation. 


The  basic  research  for  the  development  of  a  high  performance  and  high  resolution  local  polynomial  refinement 
method  in  a  steep  gradient  domain  has  been  successfully  developed  and  validated  for  combustion  simulation.  In 
a  comparative  investigation  with  a  sixth-order  compact  differencing  procedure,  the  finest  local  grid  spacing  for 
resolving  the  flame  by  the  local  polynomial  refinement  has  reached  a  value  less  7. 194x1  O'5  mm  which  is 
previously  unattainable.  The  verified  numerical  method  has  successfully  applied  to  a  counter-flowing,  stretched 
hydrogen-air  jet  flame.  The  research  results  substantiate  the  finding  for  the  criterion  of  a  minimum  flame 
thickness  for  determining  the  applicability  of  the  flamelet  concept  -  a  vortex  with  a  dimension  even  less  than 
one  millimeter  will  distort  a  thin  flame  structure  for  countering- flow  hydrogen  combustion  [3,8,9]. 

In  an  effort  to  accelerating  the  technical  transition  to  Air  Force  Research  Laboratory,  a  preliminary  investigation 
has  been  carried  out  for  combustion  within  a  hypersonic  hydrogen-air  scramjet.  This  activity  is  designed 
primarily  as  a  path-finding  effort  to  quantify  the  radiative  heat  transfer  within  the  combustion  chamber  of  a 
scramjet  with  nonequilibrium  chemical  reactions.  Additional  information  is  also  sought  in  the  scaling  of 
combustor  in  regarding  to  the  ignition  delay.  The  overall  length  of  the  simulated  scramjet  is  140  cm  with  a 
cylindrical  combustor  with  the  dimensions  of  100x17.7  cm.  The  entire  computational  domain  of  the 
axisymmetric  configuration  is  resolved  by  43,710  cells.  The  inlet  Mach  number  is  assigned  a  value  of  5.21  and 
the  fuel  is  injecting  into  the  combustor  at  a  speed  of  600  m/s  and  at  a  pressure  of  2  Atm.  The  chemical  kinetic 
model  of  the  hydrogen-air  combustion  consists  of  a  system  of  eight  elementary  chemical  reactions  for  an  un- 
premixed  supersonic  hydrogen  flame  [11].  For  the  inlet  shock  capturing,  the  Riemann  solver  is  implemented. 
The  developed  local  polynomial  refinement  therefore  is  applying  only  for  the  nonequlibrium  chemical  reacting 
flow  within  the  combustor  to  achieve  a  high  numerical  resolution. 

Figure  7  shows  the  essential  features  of  the  nonequilibrium  chemical  reactions  downstream  of  the  isolator  with 
multiple  fuel  ejectors.  The  temperature  profiles  of  the  nonequilibrium  chemical  reactions  before  and  after  the 
flame  ignited  is  presented  together  to  highlight  the  combustion  delay  and  energy  release  pattern.  The  numerical 
solutions  can  only  be  validated  by  comparing  with  theoretic  analyses  and  known  data  in  the  open  literature 
[28,29],  All  computational  results,  including  an  initiating  global  chemical  kinetic  model  reveal  the  ignition 
delay  on  the  order  of  1.32xl0‘2  seconds.  The  present  computational  result  indicates  that  the  chemical  kinetic 
model  indeed  can  adequately  describe  the  chemical  kinetic  process  and  provides  invaluable  information  for  a 
design  scaling  of  scramjet  configurations.  Furthermore,  the  complex  flame  structure  within  the  combustion 
chamber  is  captured  by  the  numerical  simulation  including  the  hot  spot  of  the  radiative  energy  exchange. 
Radiating  heat  transfer  of  the  interdisciplinary  numerical  simulation  is  described  by  both  the  diffusion  and  flat- 
layer  approximations  and  is  found  to  be  sufficient  to  model  the  additional  energy  exchange  mechanism  within 
the  combustor,  as  well  as,  to  the  interior  surface  of  the  combustor.  This  interdisciplinary  simulation  has  revealed 
the  long  overlooked  heat  protection  consideration  for  the  increasingly  higher  pressure  and  temperature 
propulsive  systems.  For  this  unique  accomplishment,  the  developed  interdisciplinary  modeling  and  simulation 
capability  is  being  transferred  to  Propulsion  and  Power  Directorate  of  AFRL  for  hypersonic  scramjet 
investigation  by  a  task  order  for  graduate  students  support. 


Acknowledgment/Disclaimer 

This  work  was  sponsored  by  the  Air  Force  Office  of  Scientific  Research,  USAF,  under  the  grant  FA9550-08-1- 
0130.  The  views  and  conclusions  contained  herein  are  those  of  the  author  and  should  not  be  interpreted  as 
necessarily  representing  the  official  policies  or  endorsements,  either  expressed  or  implied,  of  the  Air  Force 
Office  of  Scientific  Research  or  the  U.S.  Government. 

References 

1 .  Shang,  J.S.,  Computational  Fluid  Dynamics  Application  to  Aerospace  Science,  The  Aeronautical 
Journal,  Vol.  113  No.  1148,  Oct.  2009,  pp.  619-632. 

2.  Peters,  N.,  Length  Scales  in  Laminar  and  Turbulent  Flames,  in  Numerical  Approaches  to  Combustion 
Modeling,  Progress  in  Astronautics  and  Aeronautics,  Vol.  135,  1991. 

3.  Bilger,  R.W.,  The  structure  of  diffusion  Flames,  Combustion  Science  Technology,  Vol.  13,  1976,  pp. 
155-167. 

4.  Moin,  P.,  and  Hahesh  K.,  Direct  Numerical  Simulation:  a  Tool  in  Turbulent  Research,  Annual  Rev  Fluid 
Mech.,  1998,  Vol.  30,  pp.  539-578. 

5.  Pitsch,  H.,  Large-eddy  simulation  of  Turbulent  Combustion,  Annual  Rev  Fluid  Mech.,  2006,  Vol.  38  pp. 
435-482. 

6.  Grinstein,  F.  and  Kamiadakis,  G.,  Alternative  LES  and  Hybrid  RANS/LES,  J.  Fluid  Engineering,  Vol. 
124,  2002,  pp.  821-942. 

7.  Oefelein,  J.C.,  Schefer,  R.W.,  and  Barlow,  Toward  Validation  of  Large  Eddy  Simulation  for  Turbulent 
Combustion,  AIAA  J.,  Vol.  44,  No  3,  2006,  pp.  418-433. 

8.  Katta,  V.  R.,  and  Roquemore,  W.  M.,  “Simulation  of  Dynamic  Methane  Jet  Diffusion  Flames  Using 
Finite  Rate  Chemistry  Models,”  J.  AIAA  Vol. 36,  1998,  pp.2044-2054. 

9.  Katta,  V.  R.  and  Roquemore,  W.  M.,  “Calculation  of  Multidimensional  Flames  using  Large  Chemical 
Kinetics,”  AIAA  Journal,  Vol.  46,  2008,  pp.  1640-1650. 

10.  Shang,  J.S.  and  Katta,  V.R.,  Vorticity  Structure  in  Hydrogen- Air  Jet  Diffusion  Flames,  AIAA  2009- 
4036,  San  Antonio  TX,  June  2009. 

11.  Frenklach,  M.,  Wang,  H.,  Goldenberg,  M.,  Smith,  G.P.,  Golden,  D.M.,  Bowman,  C.T.,  Hanson,  R.K., 
Gardiner,  W.C.,  and  Lissianski,  V,  Gas  Research  Institute  Technical  Report  No.  GRI-95/0058,  Gas 
Research  Institute,  Chicago,  November  1,  1995. 

12.  Quealy,  A.,  Ryder,  R.,  Norris,  A.,  and  Liu,  N.-S.,  “National  Combustion  Code:  Parallel  Implementation 
and  Performance”,  NASA  TM-2000-209801,  2000. 

13.  Gautsche,  W.,  Orthogonal  polynomials,  computational  and  approximation,  Oxford  Science  Publications, 
Oxford  University  Press,  2004. 

14.  Stroud,  A.H.  and  Secrest,  D.,  Gauss  quadrature  formula,  Prentice-Hall,  Englewoodcliff,  New  Jersey, 
1966. 

15.  Korpriva,  D.,  and  Kolias,  J.,  A  conservative  staggered-grid  Chebyshev  multidomain  method  for 
compressible  flows,  J.  Comp.  Physics,  Vol.  125,  1996,  pp.  244-261. 

16.  Van  den  Abeele,  K.  Lacor,  C.,  and  Z.J.  Wang,  On  the  stability  and  accuracy  of  the  spectral  Difference 
Method,  J.  Scientific  Computing,  Vol.  37,  2008,  pp.  162- 188. 

17.  Liu,  Y.,  Vinokur,  M.,  and  Z.J.  Wang,  Spectral  difference  method  for  unstructured  grids  I:  basic 
formulation,  J.  Comp,  physics,  Vol.  216,  2006,  pp.  780-801. 

18.  Shang,  J.S.,  A  High-Resolution  Method  Using  Adaptive  Polynomials  for  Local  Refinement,  AIAA 
2010-0539,  Orlando  FL,  January  2010. 

19.  Shang,  J.S.,  Katta,  V.R.,  and  Huang,  P.G.,  Adaptive  Polynomial  local  Refinement  for  Resolving  Flame 
Structure,  AIAA  2010-4437,  Chicago  IL,  June  2010. 

20.  Roe,  P.L.,  Approximate  Riemann  solvers,  parameter  vectors,  and  difference  scheme,  J.  Comp.  Physics, 
Vol.  43,  1981,  p.351. 

21.  Shang,  J.S.,  High-Order  Compact-Difference  Schemes  for  Time-Dependent  Maxwell  Equations,  Vol. 
153,  J.  Computational  Physics  1999,  pp. 3 12-333. 


22.  Arnold,  D.N.,  An  interior  penalty  finite  element  method  with  discontinuous  elements,  SIAM  J. 
Numerical  Analysis,  Vol.  19,  1982,  pp. 742-760. 

23.  van  Raalte,  M.  and  van  Leer,  B.,  Bilinear  forms  for  the  recovery-based  discontinuous  Galerkin  method 
for  diffusion.,  Comm.  In  Comp.  Physics,  Vol.  5,  NO.  2-4,  2009,  pp.  683-693. 

24.  Huynh,  H.T.,  A  reconstruction  approach  to  high-order  schemes  including  discontinuous  Galerkin  for 
diffusion,  ALAA  2009-403,  Orlando  FL,  2009. 

25.  Gaitonde,  D.  and  Shang,  J.S.,  Optimized  Compact-Difference-Based  Finite-Volume  Schemes  for  Linear 
Wave  Phenomena,  J  Comp.  Physics  Vol.  138,  1997,  p  617. 

26.  Borges,  R.,  Carmona,  M.,  Costa,  B.,  &  Don,  W.,  An  improved  weighted  essentially  non-oscillatory 
scheme  for  hyperbolic  conservation  laws,  J.  Comp.  Physics,  Vol.  227,  2008,  p.  3191. 

27.  Lele,  S.K.,  Compact  finite  difference  schemes  with  spectral-like  resolution,  JCP.,  Vol.  103,  1992,  p.  16. 

28.  Shang,  J.S.,  Some  flow-structure  features  of  scramjet  isolator,  AIAA  2008-0772,  Reno  NV,  January 
2008. 

29.  Shang,  J.S.  and  Surzhikov,  S.T.,  Simulating  Nonequilibrium  Flow  for  Ablating  Earth  Reentry,  J.  of 
Spacecraft  and  Rockets,  Vol.  47,  No.  5,  Oct.  2010,  pp. 806-815. 

Personnel  Support  during  Duration  of  Grant 

Shang,  J.S.  Research  Professor,  Wright  State  University,  Dayton  OH 

Huang,  P.G.  Professor  and  Chairman,  Wright  State  University,  Dayton  OH 

Yan,  H.  Research  professor,  Wright  State  University,  Dayton  OH 

Surzhikov,  S.  T.  Visiting  Professor,  Russian  Academy  of  Science,  Russia 
Fabio,  R.  Visiting  Scholar,  Bologna  University,  Italy 

Andrienko,  D.A.  Graduate  student  and  visiting  Scholar 

Barnes,  C.J.  Graduate  student 

Neikamp,  M.  Student 

Publications 

1.  Shang,  J.S.,  Some  Flow-Structure  Features  of  Scramjet  Isolator,  AIAA  2008-0772,  Reno  NV,  January  2008. 

2.  Yan,  H.,  Gaitonde,  D.V.,  and  Shang,  J.S.,  Stability  Analysis  of  Thermal  Bump  in  Supersonic  Flow,  AIAA 
2008-1096,  Reno  NV,  January  2008. 

3.  Shang,  J.S.,  Huang,  P.G.,  Yan,  H.,  and  Surzhikov,  S.T.,  Electrodynamics  of  Direct  Current  Discharge,  AIAA 
2008-1101,  Reno  NV,  January  2008. 

4.  Surzhikov,  S.S.  and  Shang,  J.S.,  Kinetic  Models  Analysis  for  Super-Orbital  Aerophysics,  AIAA  2008-1278, 
Reno  NV,  January  2008. 

5.  Shang,  J.S.,  Surface  Direct  Current  Discharge  for  Hypersonic  Flow  Control,  AIAA  2008-1352,  Reno  NV., 
January  2008. 

6.  Shang,  J.S.,  Huang,  P.G.,  Yan,  H.,  Surzhikov,  S.T.,  Gaitonde,  D.V.,  Hypersonic  Flow  Control  Utilizing 
Electromagnetic -Aerodynamic  Interaction,  AIAA  2008-2606,  Dayton  OH,  April  2008. 

7.  Yan,  H.,  Gaitonde,  D.V.,  and  Shang,  J.S.,  Supersonic  Flow  Control  by  Pulsed  Thermal  Bump,  AIAA  2008- 
2607,  Dayton  OH,  April  2008. 

8.  Shang,  J.S.  and  Huang,  P.G.,  Simulating  Electrodynamic  Field  of  Dielectric  Barrier  Discharge,  AIAA  2008- 
3777,  Seattle,  WA,  June  2008. 

9.  Shang,  J.S.,  Electrostatic-Aerodynamic  Compression  in  Hypersonic  Cylindrical  Inlet,  Communication  In 
Computational  physics,  Vol.4  No.4,  2008,  pp. 838-859. 

10.  Shang,  J.S.,  Huang,  P.G.,  and  Paul,  D.B.,  Progress  in  Computational  Magneto-Fluid-Dynamics  for  Flow 
Control,  Keynote,  Fifth  International  Conference  on  Computational  Fluid  Dynamics,  July,  2008,  Seoul  Korea. 

11.  Huang,  P.G.  and  Shang,  J.S.,  High-Resolution  Simulation  of  Direct  Current  Discharge  Coupled  with  a 
Magnetic  Field,  Proceeding  of  XVII  International  Conf.  on  Gasdischarge  and  their  Application,  Cardiff,  UK, 
Sept.  2008,  pp.441-444. 

12.  Shang,  J.S.,  Kimmel,  R.L.,  Menart,  J.,  Surzhikov,  S.T.,  Hypersonic  Flow  Control  Using  Surface  Plasma 
Discharge,  J.  Propulsion  and  Power,  Vol.  24,  No.  5,  2008,  pp.923-934. 


13.  Shang,  J.S.,  Surface  Direct  Current  Discharge  for  Hypersonic  Flow  Control,  J.  of  Spacecraft  and  Rockets, 
Vol.  45,  No.  6,  2008,  pp  1213-1222. 

14.  Shang,  J.S.,  Surzhikov,  S.T.,  and  Yan,  H.,  Simulate  Hypersonic  Nonequlibrium  Flow  Using  Kinetic  Models, 
AIAA  2009-0386,  Orlando  FL.  January  2009. 

15. Shang,  J.S.,  Huang,  P.G.,  Yan,  H.,  Surzhikov,  S.T.,  Computational  Simulation  of  Direct  Current  Discharge, 

J.  Applied  Physics,  Vol.  105,  2009,  pp.  023303-1-14. 

16.  Shang,  J.S.  and  Katta,  V.R.,  Vorticity  Structure  in  Hydrogen- Air  Jet  Diffusion  Flames,  AIAA  2009-4036, 
San  Antonio  TX,  June  2009. 

17.  Shang,  J.S.,  Analysis  of  Hypersonic  Nonequilibrium  Flow  over  Ablating  Surface,  AIAA  2009-4051,  San 
Antonio  TX,  June  2009. 

18.  Shang,  J.S.,  Huang,  P.G.,  and  S.  Stanfield,  Electrodynamic  Field  of  Dielectric  Barrier  Discharge,  17th 
International  Conference  on  MHD  Energy  Conversion,  Kanagawa,  Japan,  14-17  July  2009. 

19.  Shang,  J.S.,  Computational  Fluid  Dynamics  Application  to  Aerospace  Science,  The  Aeronautical  Journal, 
Vol.  113  No.  1 148,  Oct.  2009,  pp.  619-632. 

20.  Shang,  J.S.,  A  High-Resolution  Method  Using  Adaptive  Polynomials  for  Local  Refinement,  AIAA  2010- 
0539,  Orlando  FL,  January  2010. 

21.  Huang,  P.G.,  Shang,  J.S.,  and  Stanfield,  S.A.  Periodic  Electrodynamic  Field  of  Dielectric  Barrier  Discharge, 
AIAA  2010-0961,  Orlando  FL,  January  2010. 

22.  Shang,  J.S.  and  Surzhikov,  S.T.,  Simulating  Ablation  Phenomenon  for  Earth  Reentry,  AIAA  2010-0983, 
Orlando  FL,  January  2010. 

23.  Shang,  J.S.  and  Huang,  P.G.,  Modeling  AC  Dielectric  Barrier  Discharge,  J.  Applied  Physics,  Vol.  107, 
No.ll.pp  113302-1-7.  2010. 

24.  Shang,  J.S.,  Katta,  V.R.,  and  Huang,  P.G.,  Adaptive  Polynomial  local  Refinement  for  Resolving  Flame 
Structure,  AIAA  2010-4437,  Chicago  IL,  June  2010. 

25.  Surzhikov,  S.T.  and  Shang,  J.S.,  Coupled  Radiation-Gasdynamics  Model  for  Stardust  Aerothermodynamics 
Data,  AIAA  2010-4521,  Chicago  IL,  June  2010. 

26.  Shang,  J.S.  and  Surzhikov,  S.T.,  Study  Physical-Chemical  Kinetic  Models  for  Simulating  Stardust  Earth 
Reentry,  AIAA  2010-4455,  Chicago  IL,  June  2010. 

27.  Shang,  J.S.  and  Huang,  P.G.,  A  local  resolution  refinement  algorithm  using  Gauss-Lobatto  quadrature,  the 
sixth  international  conference  on  Computational  Fluid  Dynamics,  St.  Petersburg,  Russia,  July  12-16,  2010. 

28. Shang,  J.S.  and  Surzhikov,  S.T.,  Simulating  Nonequilibrium  Flow  for  Ablating  Earth  Reentry,  J.  of 
Spacecraft  and  Rockets,  Vol.  47,  No.  5,  Oct.  2010,  pp. 806-815. 

29.  Shang,  J.S.,  Analytical  Foundation  of  Hypersonic  Flow,  Encyclopedia  of  Aerospace  Engineering,  John  Wily 
&  Sons,  Ltd,  England,  2010 

30.  Barnes,  C.J.,  Huang,  G.P.,  and  Shang,  J.S.,  A  high  resolution  spectral  method  for  the  compressible  Navier- 
Stokes  Equations,  AIAA  2011-0049,  Orlando  FL,  January  2011. 

31.  Surzhikov,  S.  and  Shang,  J.  Radiative  and  convective  heating  of  Orion  space  vehicles  at  Earth  orbital 
entries,  AIAA  2011-0251,  Orlando,  FL,  January  2011. 

32.  Shang,  J.S.  and  Surzhikov,  S.T.,  Nonequilibrium  radiative  heat  transfer  in  hypersonic  flow,  AIAA  2011- 
0762,  Orlando,  FL,  January  2011. 

33.  Fabio,  Roveda,  Shang,  J.S.,  and  Huang,  P.G.  On  the  electrodynamic  force  of  dielectric  barrier  discharge, 
AIAA  201 1-0905,  Orlando,  FL,  January  2011. 

34  Huang,  P.G.,  Shang,  J.S.,  and  Stanfield,  S.A.,  Periodic  Electrodynamic  Field  of  Dielectric  Barrier 
Discharge,  AIAA  J.  Vol.  49,  No.l.,  January  2011,  pp.  119-127. 

35.  Shang,  J.S.  and  Surzhikov,  S.T.,  Nonequilibrium  radiative  heat  transfer  in  hypersonic  flow,  AIAA  2011- 
2258,  San  Francisco,  CA,  April,  2011. 

36.  Shang,  J.S.  and  Surzhikov,  S.T.,  A  prospective  of  physical  modeling  for  nonequilibrium  hypersonic 
simulation,  7th  European  Symposium  on  Aerothermodynamics  for  space  Vehicles,  Brugge  Belgium,  May  2011. 

37.  Surzhikov,  S.T.  and  Shang,  J.S.,  Radiative  aerothermodynamics  of  entry  probes  in  Martian  and  Earth 
atmosphere,  7th  European  Symposium  on  Aerothermodynamics  for  space  Vehicles,  Brugge  Belgium,  May 
2011. 

38.  Shang,  J.S.,  Roveda,  F.,  and  Huang,  P.G.,  Electrodynamic  force  of  dielectric  barrier  discharge,  J.  Applied 
Physics,  Vol.  109,  No.l  1,  pp.  113301-8,2011 


39.  Shang,  J.S.  and  Surzhikov,  S.T.,  Simulating  Stardust  Earth  Reentry  with  Radiation  Heat  Transfer,  J.  of 
Space  craft  and  Rockets,  Vol.  28,  No.  3,  201 1,  pp.385-396. 

40.  Shang,  J.S.,  Katta,  V.R.,  and  Huang,  P.G.,  Adaptive  multi-grid  blocks  polynomials  refinement  for  flame 
structure  resolution,  submitted  to  Communications  in  Computational  Physics. 

Honors  &  Awards  Received 

USAF  Basic  Research  Award,  1986 
Fellow  ofAIAA,  1993 

Outstanding  and  Exception  Civilian  Service  Awards,  2001 
AIAA  Plasmadynamics  and  Laser  Award,  2004 

AFRL  Point  of  Contact 
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Dr.  Richard  Rivir,  AFRL/RZ  WPAFB,  OH  937-255-2744.  Sponsor  for  the  FA8650-08-D-2806,  Task  Order  4, 
A  path- finding  study  for  radiative  exchange  in  combustion  chamber,  interacted  monthly. 

Mr.  Charles  Brinks,  USAF  X51A  Program  Manager,  AFRL/RZAT,  WPAFB,  937-255-7611 
Dr.  Roger  Kimmel,  AFRL/RBAA  WPAFB,  OH  937-255-8295,  interacted  weekly. 

Mr.  Michael  Zeigler,  AFRL/RBAI  WPAFB,  OH  937-656-6307,  interacted  and  met  weekly. 

Dr.  Melvyn  Roquemore,  AFRL/RZ  WPAFB,  OH  937-255-68 1 3 
Dr.  Biswa  Ganguly,  AFRL/RZ  WPAFB,  OH  937-255-6782 

Transitions 

The  present  basic  research  was  initiated  to  introduce  innovative  numerical  algorithms  and  to  develop  high- 
resolution,  interdisciplinary  simulation  capability  for  combustion.  Based  on  this  long-term  objective,  the 
technical  transitions  are  focused  on  the  dissimilation  of  research  results  first  to  the  research  scientists  of  AFRL 
and  then  to  national/intemational  conferences  of  professional  societies.  All  the  technical  transition  processes  are 
consistently  followed  up  through  personal  interactions. 

The  most  effective  technical  assist/transition  activities  can  be  summarized  as  (1)  on  September  2011,  Maj.  Gen. 
W.N.  McCasland,  Commander  of  Air  Force  Research  Laboratory,  invited  the  principal  investigator  to  serve  on 
the  External  Independent  Review  Team  (EIRT)  for  performance  views  of  the  USAF  X51A  Scramjet  system.  (2) 
The  AFOSR  sponsored  local  grid  refinement  computational  simulation  for  nonequilibrium  combustion  is 
expanded  its  applicability  from  combustion  to  include  radiative  heat  transfer  for  advanced  propulsion  systems. 
The  basic  research  team  received  a  task  order  (FA8650-08-D-2806,  Task  Order  0004)  for  graduate  student 
support  from  the  Propulsion  Directorate  of  the  Air  Force  Research  Laboratory  to  create  a  path-finding 
procedure  to  study  the  radiative  exchange  within  combustor. 

The  additional  technology  assist  effort  to  USAF  includes  the  transition  of  the  concept  and  potential  practical 
applications  of  the  Gaussian  Quadrature  to  premixed-  and  non-premixed  turbulent  hydrogen  combustion.  A 
scientific  interaction  has  been  maintained  with  Dr.  W.  Melvyn  Roquemore  (AFRL/RZ  937-255-6813).  The 
continuous  exchange  and  collaboration  has  also  been  sustained  with  the  combustion  research  group  of  the 
Propulsion  Directorate.  This  activity  has  proved  to  be  extremely  mutual  benefit  in  advancing  basic 
understanding  of  extremely  complex  combustion  phenomena  for  USAF’s  interest. 

The  new  model  and  computational  simulation  for  combustion  have  also  been  conveyed  to  Dr.  Biswa  Ganguly;  a 
mutual  support  agreement  for  plasma  assisted  ignition  and  enhanced  combustion  stability  has  been  reached  in 
January  2010.  The  present  project  will  provide  a  direct  technical  support  to  Center  of  Advanced  Power  and 
Energy  Conversion  of  the  joint  Wright  State  University  and  Air  Force  Research  Laboratory  (AFRL/RZPE  937- 
255-6782). 

New  Discoveries 

None. 


